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^ Abstract 

We analyze parity-violating nucleon-nucleon scattering at energies using the effective 

^ field theory appropriate for this regime. The minimal Lagrangian for short-range parity-violating 

\^ NN interactions is written in an operator basis that encodes the five partial-wave transitions that 

dominate at these energies. We calculate the leading-order relationships between parity- violating 
I— I NN asymmetries and the coefficients in the Lagrangian and also discuss the size of sub-leading 

corrections. We conclude with a discussion of further observables needed to completely determine 

the leading-order Lagrangian. 
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I. INTRODUCTION 



The existence of parity violation in nuclear forces is a manifestation of the presence of 
weak interactions between the quarks in the nucleon. In this paper we discuss the most 
basic observables that display this phenomenon: A^A^ scattering asymmetries that would 
be zero were parity conserved in the A^A^ interaction. We do this using an effective field 
theory (EFT) that is based on the existence of large l/'n^vr) scattering lengths in the A^A^ 
system [Il[2l[3llll|5l|6l[71[Hl[9l HD] • The scattering experiments we will study allow us to use 
an EFT that contains only nucleon degrees of freedom, and treats all interactions between 
those degrees of freedom as being short-ranged. We choose to consider only the scattering 
length anomalously large; higher order corrections are obtained as a systematic expansion 
in powers of r/a, where a is the two-body scattering length and r the effective range in the 
corresponding partial wave. 

We find it convenient to examine such corrections to parity- violating A^A^ observables us- 
ing a dynamical dibaryon field [HI [121 US] • This treatment captures the dominant dynamics 
in the situation where both the scattering length and effective range are unnaturally large. 
It corresponds to resumming an infinite subset of terms that are higher order in the power 
counting we use, where r ~ 1/^,7^ ^ a ~ Therefore, expanding the results obtained 
from the dibaryon formalism in powers of r/a reproduces the case of interest to us here, and 
so we employ dibaryon fields as a calculational tool. 

There are a number of theoretical treatments based on hadronic degrees of freedom that 
have been used to study parity violation (PV). Pioneering theoretical studies on PV in 
nuclear forces were carried out by Danilov [13] and Desplanques and Missimer [T^]. For 
most of the last thirty years, the framework of single-meson exchange, most commonly 
using the taxonomy developed by Desplanques, Donoghue, and Holstein (DDH) [16], has 
been the one used to interpret and motivate experiments [T71 [T8l [T9] . But the lack of a 
concordance region in the space of DDH PV parameters (see, e.g., the plot in Ref. [18]) may 
be related to the model assumptions — such as the mediation of the PV A^A^ interaction by 
vector mesons — in that approach. The consequent desire by the community to "...recast the 
DDH language.... in terms of... effective short-range parity- violating N-N interactions" [20] 
has motivated EFT treatments of the PV A^A^ force. 

Recently the version of chiral perturbation theory appropriate for few-nucleon systems [H, 
EH |22l |23l EH |25] (xET hereafter) has been used to derive the long-range (r ~ l/m^r) part of 
the PV A^A^ force, and to classify the short- distance operators that appear in that force [26] . 
(See also the original xPT analysis of PV pion- nucleon operators in Ref. [27|). This has 
the advantage that — up to a given order in xET — one can guarantee that a complete set 
of parity-violating operators has been considered. The potential of Ref. [26] is formulated 
in terms of the appropriate degrees of freedom for momenta of order m.,^: nucleons and 
pions.^ The use of heavier mesons to encode A^A^ interactions of different t-channel quantum 
numbers is not necessary at energies below ~ 200 MeV. By fitting the constants that encode 
the short-distance A^A^ interaction to data, the xET treatment of PV in A^A^ scattering 
avoids any assumption about what dynamics is at work for r <C The consequences of 

xET for PV A^A^ scattering, as well as other PV few-nucleon-system observables, have been 
computed in Refs. [30l |3Tl |32]. However, the PV operators and strong-nuclear- force wave 
functions employed in these works were not consistent, since phenomenological models were 



^ For extensions to include the A isobar see Refs. [28l [29]. 



2 



used for the latter, but a xET for the former. Refs. [SOj |32] also include what are referred to 
as "pionless-theory" results, and there the mismatch between operators and wave functions is 
a serious problem. The use of AV18 wave functions for one, and a pionless EFT for the other, 
involves a mismatch of roughly an order of magnitude in the resolution (~ renormalization) 
scales of these two different calculational ingredients. 

The pionless theory, EFT{jt), is relevant to studies of parity violation because many 
of the existing and planned experiments [33l IMl EHl EEl |37l [Ml EH], HO] take place at 
energies at or below lO's of MeV. In this region the pion-exchange nature of the nuclear 
force is not resolved. For energies E < rn^/M the EFT in which interactions are encoded 
as A^A^ contact operators contains all the relevant degrees of freedom. The convergence 
of this EFT for two-nucleon-system observables is well demonstrated O |10l HH |32] in the 
low-energy regime. In contrast, there have been significant questions raised recently about 
the appropriate power counting for short-distance operators in the NN system in the xET 
where pions are explicit degrees of freedom [5], HSl HH HHl HGJ HT] , as well as about whether 
the Delta(1232) needs to be included as an explicit degree of freedom in order to guarantee 
reasonable convergence [211 HHl 09] • 

In EFT(/'), parity violation in A^A^ scattering is described by contact operators that 
have a lowest possible dimension of seven. As Girlanda has recently shown [50] there are 
five such independent operators.^ These five leading-order EFT contact operators contain 
the same physics as the five Danilov amplitudes [T^ that encode the mixing between S 
and P-waves that becomes possible in the presence of parity violation. In Section [IT] we 
rewrite the Lagrangian of Ref. [HD] in terms of five operators that each mediate a specific S-P 
transition. This makes the computation of longitudinal asymmetries in A^A^ scattering (given 



in Sec. Ill) straightforward. We present analytic results for the longitudinal asymmetries 



in nn, pp, and np scattering. We use a dynamical dibaryon field to obtain a portion of the 
higher-order corrections in r/a and rp in the EFT(^) expansion for the strong rescattering. 
Subsection IIIB computes the Coulomb effects that are present in pp scattering at low 
energy and so modify the expression in the pp case. We extract the scale dependence of the 
result obtained when assuming that the r dependence is higher order, and then discuss the 
r dependent corrections. In Sec. llVlwe compare our results with the two existing pieces of 



experimental data [55| IM] that are within the range of validity of this EFT. We close in 
Sec. |V]with a summary and a discussion of further A^A^ system experiments that could pin 
down the LO PV EFT(;f) Lagrangian. Details of the conversion of one set of operators to 
another are given in an Appendix. 



II. LAGRANGIANS 

At low enough energies, the details of the gauge boson ((yf, W^, Z) exchange between 
interacting quarks in the two-nucleon system are not experimentally accessible. Instead, 
the system can be described by using nucleon interpolating fields and treating both strong 
and weak interactions as contact interactions. To a given order there are a finite number of 
independent operators that describe these interactions. 



^ While Ref. [26^ lists ten operators, it is pointed out that only five different combinations are relevant at 
low energies. 
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The Lagrangian for the parity-conserving two-nucleon sector may be written as [T] 



^ = ^^(^^0 + - l^siN^^y - ICriN^aN)' + ..., (1) 

where the elhpsis stands for terms that contain more derivatives, and the nucleon field N 
carries both isospin and spin indices. The cxj are the SU(2) Pauli matrices in spin space 
and the will be the SU(2) Pauli matrices in isospin space. In Eq. ([T]) the A^A^ effective 
ranges are assumed to be "natural" compared to the expected scale of 1/171,^ while the A^A^ 
scattering lengths are large in the same units. This facilitates an expansion of the A^A^ 
amplitude in powers of the small parameter Q, where Q ~ 1/a ~ p (with p the A^A^ relative 
momentum) |2l |3l E] . 

The same physics can be expressed using an operator basis that makes the incoming and 
outgoing partial waves explicit [5T| [52]. 

C =N^{zdo + - ^Cr«)(Ar^r2r,a2Ar)t(Ar^r2r,a2Ar) 

- lcj,'''\N^T2a2aiN)\N''T2a2aiN) + ..., (2) 

o 

where now Cq = Cs — SCt and Cg = Cs + Ct- The operator between nucleons is simply 
the projector onto the relevant partial wave, with the normalized projectors being [SU [52] 

PaCSo) = -^T2Taa2 ] Pi^Si) = \.T2(J2(Ji ■ 



A convenient form of this Lagrangian is provided by use of dibaryon fields pTl 
This form is equivalent to Eqs. ^ and ([T]) at leading order, but resums all the higher-order 
corrections in A^A^ scattering that are proportional to the effective range. It would therefore 
give the exact A^A^ amplitude in both the ^Sq and channel were the shape parameter 
and all higher-order terms in the effective-range expansion zero. Dibaryon fields, Sq and tj, 
respectively, for the ^Sq and ^5*1 states, are included in the Lagrangian [T3l [53] : 

>C = Art + j iv - 4 (^^Oo + ^ - A(35,) j U - 9es,) [tlN^P^{%)N + h.c. 

-4 (^^do + ^ - A(.5„) j - g^rso) [siN^PaCSo)N + h.c] , (3) 

with a an isospin and i a spin index. If PDS (power divergence subtraction) [3l [1] is used 
to compute loops then for both channels we have (in an obvious notation with channel 
subscripts suppressed) : 

The leading-order PV two-nucleon Lagrangian can also be expressed in a variety of bases. 
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Paralleling the one used in Eq. ([T|, Ref. |50] writes: 

-^ie.,fciVV,iVV,(iVtafciV) 

- G2e^jk [NW^aiNVjiN^akN] + N^aiNVjiNW^akN)] 

- gaabe^JkN^TaCTiNVJ{N^n(TkN) 

+ GeeabsViN^TaN) ■ N^naN^ , (5) 

<->■ ^ — 
where a V & = ciVb — a'^h and 

1 = 




Note that we have renamed the coefficients as compared to Ref. [50], in order to avoid 
confusion with the parity-conserving Lagrangian. In doing so we absorbed into the ^s the 
overall normalization factor of that multiplies the Cj's and Ci 's in Ref. [50]. In 

EFT(^) the coefficients are typically dependent on the renormalization point, /i, used in the 
evaluation of loop diagrams. Using the partial- wave basis, as in Eq. ([2]), we have: 



PW 
PV 



C(2:oT°^ {N^(y2r2rN)^ (iV^^2 o-tV t^tN'^ + 
C^iZT {N^^2T2r''N)^ {N^a, a- V r^r^iv) + 
<^(a1=2T°^ [^^^2r2r^^)^ [^^^2 o-iV r^r'N) + 

The two Lagrangians in Eq. Q and Q give the same results for physical observables if 
the low-energy constants obey the relationships (see Appendix): 









'-'(A7= 


=0) 




'-'(A7= 


= 1) 




'-'(A7= 


-^Po) 
=2) 


= -2^5 , 




-3Pi) 


= i^6 • 



(7) 



Note that there are five independent coefficients at this order — as explained in Refs. [261 [50l 
51]. They dictate the only possible nucleon-nucleon scattering observables at low enough 
(i.e., non-dynamical pion) energies. From the partial- wave point of view, only the coefficients 
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£{^Si ipi)^ ^{Ai=o)'°\ ^(Ai=2^°^ y ^^^^^ involved in parity- violating neutron-proton 

observables; while ^{ai=i)'°^ ^ ^^'^ ^{ai=2)"^ involved in parity- violating neutron- 

neutron (or proton-proton) observables. 

In the DDH approach, all but the ^Si — ^Pi operator are considered only in terms of 
vector-meson exchange. At low energies, the vector mesons are not dynamical, so the DDH 
description can be considered a way to "encode" the processes so that calculations and exper- 
iments can be compared, so long as the vector-meson interpretation is not taken literally. In 
particular, at very low energies only five independent parameters are relevant for the physics 
of parity violating NN scattering — they are sufficient to encode all leading-order phenom- 
ena. The EFT parameterization presented here provides a model-independent language in 
which to compare experiments. 



III. RESULTS 

A. Calculation of longitudinal analyzing power 

Parity violation in the NN interaction leads to mixing between odd and even partial 
waves. To obtain the leading effects of this mixing in EFT{-/i) it is sufficient to calculate the 
amplitude that mediates S-wave to P-wave transitions. The mixing of higher partial waves 
is suppressed by additional powers of the small parameter Q. 

The leading diagrams contributing to this parity-violating A^A^ scattering amplitude are 
shown in Fig. [1} Note that only the S'-wave side receives an enhancement from the strong 
S-wave bubble sum. Diagrams with strong rescattering on the P-wave side are higher order. 
We evaluate the diagrams shown in Fig. [l] using the PDS [21 S] renormalization scheme to 



calculate the loops. Keeping in mind the issue of higher-order corrections (see Subsec. IIIC) 
we employ the dibaryon Lagrangian ^ to compute the strong rescattering. Our result for 
the scattering amplitude is: 



Tr4, = ±4p An/pp ( ^ - V^V^ -^^ (-^^- + ^p] , (8) 





FIG. 1: Diagrams contributing to parity- violating NN scattering. The shaded blob is the leading 
order parity-conserving amplitude, while the square denotes the vertex from the parity-violating 
Lagrangian. The nucleons are in a P-wave on the left of the parity-violating vertex, and in an 
S-wave on the right. 
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and 



PV 

np 



±4p 



1 



1 



+A. 



np 



1 



r ^°p'^ + ip 



1 1 3cj 2 



-1 



(9) 



where —iT is the sum of diagrams in Fig. [T| p = |p| (see Fig. [2]) and the upper (lower) sign 
is for a beam of positive (negative) hehcity. The (strong) parameters a^^^^-' and r^^^^^J 
are, respectively, the scattering length and effective range of a particular partial wave. The 
weak A^A^ interaction parameters are collected in amplitudes A, which are given by: 



Ann = Gl + Gl — 2(^2 + G5 



4 C 



'-'(A/=l) ^(AI=2) 



A., 



VP 



,(iSo-3Po) 
^(A7=0) 

G1 + G1 + 2{G2 - G5) 

^ ' '-'(A/=0) '-'(A/=l) ^{AI=2) 



Alf; = Gi + Gi + 4G5 



4 C 



,(i5o-3Po) 



(10) 

(11) 

(12) 
(13) 
(14) 
(15) 
(16) 
(17) 



in the notation of the Lagrangians ([5]) and ([6]), respectively. While the expressions 
appear to have an explicit subtraction point (/i) dependence, as physical observables each 
must be /i independent. This dictates the scaling of the Ann with respect to fi. 

The leading-order [0((5°)] amplitudes are obtained by setting the effective ranges r ^° 
and r^^^ equal to zero in Eqs. (jsj) and (9|. This yields, for example. 



A. 



np 



^^(i5o-3Po) 
'(A7=0) '^'-'(A/=2) 

Gl - Gl - 2^6 

4 (c^''S'-'^''> - 2C('^i~'^i) 



nn V 



0) = ±4p 



Ann 47r 



(18) 



a result that already appeared in Section 4 of Ref. |26], but here we have also provided the 
relationship of the parity-violating amplitude Ann to the coefficients in the Lagrangian(s) 

(Eqs. (10) and (11)). The ratio of Ann to Cq ° must be independent of /i. Since, in PDS, 



a 



{'So) 



Air 



M4r 

a So 



we have 



1 OA nr 

Ann dfl 



1 

a'^0 



/i 



(19) 



(20) 



In NN scattering at these energies the weak interaction is about 10 times the strong 
interaction. Therefore feasible experiments involve observables that vanish under strong 
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interactions. Relevant A^A^ measurements have focused on longitudinal asymmetries in iV + 
N scattering. Here, the interference terms between the strong and weak operators change 
sign when the longitudinal polarization of the incoming nucleon changes sign. The strong- 
interaction scattering is unaffected by a change in polarization, so an asymmetry is formed 
when the differential cross sections of the two different polarization states are subtracted. 

From the scattering amplitude calculated above we can determine the longitudinal asym- 
metry: 

= 

a+ + cr_ 

where a± is the total scattering cross section of a nucleon with helicity ± on an unpolarized 
nucleon target — unless integration over a restricted angular range (e.g. in pp scattering) is 
indicated. Neglecting, for the moment, the Coulomb interaction in the pp case we find 



2M 



IT 



-pAr. 



1 



1 



T^^"p'^ 



(22) 



and 



2M 



TT 



pAp 



1 1 



r - /i I , 



(23) 



AT 



2M 



IT 



-p 



dQ 

dn 



■A. 



'So 



dn 



A 



dn 



dn 



1 1 



1 1 3'?, 2 

3sr - o"" P 



(24) 



The differential cross sections and ^jj^ 
conserving Lagrangian (see Eqs. Q and (|2|): 

da _ fl 1 



only contain contributions from the parity- 



n -1 



2 \ 2 

rp ] + p 



(25) 



We have again suppressed the channels' superscripts. 

Upon setting r "^^ = and using Eq. (19), Eqs. (22) and (23) recapture the form derived 
in Refs. [261 [M]: 



.pp/nn 



A. 



■pp/nn 



'-'0 



This, and the more complex formula for A^'- 



(26) 



A7 



•ip 



dn 



dn 



"^np 



+ 



dn 



dn 



+ 3 



da' 



•^np 



dn 



(27) 



are the LO predictions of EFT{-/i) for these asymmetries. The only unknown quantities 
in these predictions are the coefficients of the parity- violating Lagrangian. Eqs. (22)-(24) 



relate these coefficients (Eqs. (10)-(17)) to observable asymmetries. From these expressions 
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we see that a measurement of all three analyzing powers as a function of energy could pin 
down four different combinations of coefficients, since the two pre-factors in Eq. (24) have 
distinct energy dependence — even if r^^° = r^'^^ = 0. However, the nn experiment is not 
feasible in the foreseeable future, and so alternative strategies to access the combination Ann 
are probably necessary. 



B. Coulomb corrections for pp 



The result in Eq. (26 ) ignores the Coulomb interaction. Coulomb photons can be included 
in EFT(^), and the computation of Coulomb scattering was carried out to leading order for 
S-wave NN scattering in Ref. [SS, EB] . 

In the parity- violating case the computation proceeds as in Fig. [T} except now Coulomb 
photons must be added to the initial, final, and all intermediate states. Since the initial- 
state and final-state Coulomb scattering factorizes this yields the final result, quoted in 
Ref. 1261 [S]: 



A 



NC\ 



0) = ±4p-f C,^exp(^(ao(r/) + aM))^ 

^0 



1 



Up) 



(28) 



where the purely Coulombic part of the scattering amplitude has been separated off the 
total amplitude [57] . 

T = Tnc + TcouI. (29) 



is the Sommerfeld factor: 



^2 



(30) 



with the Coulomb parameter r] = -i^^ and (J[{r\) = argr(Z + 1 + irj), where T is the Euler 
gamma function. Jo(p) is the Coulomb-modified bubble: 



An 



H{r]) - In 



aM 



An 



once divergences in D = A and D = 3 have been subtracted, 

H{r]) = ip{ir]) + 7^ - log(i?7). 



(31) 



(32) 



with ip the derivative of the Euler Gamma function, and Ce = 0.5772(...) is Euler's constant. 
(See also Refs. |MllSn].) 

Experimental asymmetries are typically measured over a finite angular range. This is a 
particularly important detail in pp scattering, due to the infinite Coulomb cross section in 
the forward direction. Implementing the integrals over a finite range 61 < 6 < 62 we have: 



de sine 2Re[(T^^)(T^c + TcouiV] 



Af 



Jg^ dO sine \Tnc + Tcoui\ 



2 Re 



rpPV 

-'- NC 



T, 



NC 



COS 6*1 — cos 6*2 Jq^ 



02 



dO sin 9 



T, 



Coal 



NC 



(33) 
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where we have used the fact that Tnc,pv and T^c are angle independent at this order, and 
have neglected |Tco„;p. Ref. [55] finds that 



T, 



Q2 g2jcro{»?) 



NC 



yielding 



Af 



■ip 



A, 



VP 



Re 



^j[(Ti(-^)-(To(»y)] 



COS 6*1 — cos 6*2 



02 



dO sin 9 



T, 



Coul 



NC 



(34) 



The factor in square brackets contains the Coulomb corrections to the result of the previous 
section. It is a function of ?7, the scattering length a, and the angular range being examined. 
For the experiments of interest here we have r] For small r], 



Coul 



and 



T, 



1 — cos 6 

4:71 1 



NC 



M 



with the strong pp scattering length, as{fi), defined by: 

47r 1 



-/i. 



(35) 
(36) 

(37) 



A reliable extraction of as{fi) from pp data appears to require a computation to several 



orders in EFT(/') ^60j. Instead, for comparison to experiments in Sec. IV, we use isospin 



symmetry, and take for as{jj,) the 'recommended' central value of the nn scattering length, 
Onn = —18.59 fm [61]. We obtain 



ip 



'-'0 



Re 



ip 



■pp 



l + V 



l + V 



as{n)p 



+ i 



cos 6*1 — cos 6*^ 



In 



cos 6*1 



cos 6*2 



0's{iAp J cos Q\ — cos 6*- 



In 



cos 6^1 



cos 6*5 



(38) 
(39) 



so long as forward angles are avoided. 

Even for pp experiments at Tiab = 0.5 MeV, we have ?7 ~ 0.22, so ri should be a good 
expansion parameter. Since parity-violating asymmetries grow as \fT\^ the extant measure- 
ments of Ai^ were conducted at energies significantly higher than this, so in practice ignoring 
Coulomb (as was done in the pioneering study of Ref. [IS]), or expanding in powers of r], is 
a good approximation. 

This suggests using a different expansion where effects proportional to Ma/p are treated 
perturbatively. However, numerically Ma ~ 1/a in the ^5*0 channel. If Ma/p is treated as 
a small parameter, then \/[ap) should really also be treated as a perturbation. This results 
in a theory set up as an expansion around the unitary (|a| — > oo) limit [62j. Attempts to 
treat Coulomb interactions in perturbation theory and retain the 1 / (ap) corrections in the 
unitary limit to all orders requires care since the divergences that are present in the Coulomb 
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bubble must still be absorbed [551 [63]. While the Coulomb contributions to the final result, 
Eq.(39), are small for all existing and proposed experiments, here we have retained them 
to all orders in the intermediate steps of the calculation of the pp scattering amplitude, 
and only performed the expansion in powers of ry when that amplitude is inserted in the 
expression for the asymmetry. 



C. Corrections proportional to the effective range 



Here we discuss corrections to the leading-order result of Eq. (26). Since we only examine 



the form of the NLO correction, and do not compare to experimental data, we consider nn 
scattering. The arguments are similar for pp scattering and np scattering. Only operators 
with the same space-spin structure as the leading-order ones of Eq. ^ are necessary for this 
analysis. Other space-spin structures, e.g, mixing between P- and D-waves, have the same 
number of derivatives as the operators we will consider in this section, but the resulting 
amplitudes are not enhanced by the strong S-wave rescattering. (E.g, effects of P-D mixing 
do not enter until 0{Q^): three orders beyond leading.) 

The result in Eq. ( [22| is actually somewhat deceptive. The use of the dibaryon formalism 
in the strong Lagrangian seems to imply that the physics of the effective range has been 
included to all orders in Eq. (22). This is not the case because the dibaryon formalism 
was not used in the weak Lagrangian. Either scaling for the effective range (as or Q^) 
will lead to consistent results — but only if the choice is used uniformly in all aspects of the 
calculation. 



In particular, demanding that 



Q^nn 



(40) 



implies that the Ann of Eqs. (10) and (11) becomes energy dependent: 



A 



fx 



(41) 



where we have dropped the partial-wave specification on a and r. To obtain consistent results 
for the case considered here, where the effective range is natural (~ Q°), we must expand 



Eq. (41) in powers of r. This allows us to estimate the impact of corrections proportional 



to r in the weak-interaction piece of the Lagrangian. It yields: 



A 



1 



1 + 



1 rp^ 



Writing this as 



A 



aLO I -2 AISLU I 



NLO 



+ ■ 



(42) 



(43) 



makes it clear that there must be corrections to the leading-order weak-interaction La- 
grangian that have the same space-spin structure, but are proportional to the square of the 
momentum (equivalently, the energy) of the NN collision. Each term will be accompanied 
by its own low-energy constant. These corrections to the Ann of Eqs. (10) and (11) are sup- 
pressed by one power of the small parameter Q. They are missing from the result (Sj), which 
includes only the leading-order part of Ann, and the effect of strong rescattering. Neglecting 
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NLO contributions to the weak Lagrangian results in an inconsistent calculation as soon as 
r ^ 0. 



While A^^ recaptures the scaling of Cq (Eq. (19)), as expected, A^^'^ runs like C2, the 



NLO strong coefficient (see, for example, Eq. (2.26) of Ref. |1]): 

aNLO 1 



'■nn 



(44) 



The necessity for the weak Lagrangian to have an 0{Q) piece with coefficients scaling ac 



cording to Eq. (44) can also be derived by considering the /x-invariance of the NLO amplitude 
for PV A^A^ scattering in a strictly perturbative calculation in powers of Q. Conversely, were 
we to use a weak Lagrangian expressed using dibaryon fields for the S'-channels the scaling 



(41) would emerge automatically. In either case, in order to maintain /i-independent results 
the counting of r must remain consistent between Cweak and Cstrong- Both Lagrangians 
contain higher-order terms that are proportional to r, and the scaling of these contributions 
with fi is correlated. If effects proportional to r are resummed using a dibaryon formalism 
consistently in both weak and strong Lagrangians, naive dimensional analysis suggests that 
additional corrections in Cpv, which are related to additional parameters, are suppressed 
by two powers of Q. 



IV. COMPARISON WITH EXPERIMENT 

In order to completely specify the leading-order PV NN Lagrangian in this EFT the 
coefficients of the five dimension-7 operators must be determined. The only way to do this 
in a model-independent fashion is to fit them to experiment. If the experiments are at low 
enough energy the corrections from higher-dimensional operators that encode other partial- 
wave transitions, as well as energy-dependence of the S-P transitions, will presumably be 
small. In practice a higher-order analysis, together with a variety of different measurements, 
will have to be employed to see if the EFT is complete and consistent. 

Here we pursue only a leading-order analysis of the two most recent low-energy measure- 
ments of the longitudinal asymmetry in pp scattering. These yielded 



Af{E = 13.6 MeV) = (-0.93 ± 0.21) x 10"^ (45) 

and 



A^{E = 45 MeV) = (-1.50 ± 0.22) x 10"^ (46) 
in the angular range 23° < 9iab < 52°. 



Using (39) together with the lower-energy number (45) yields: 



Appifi = m^) = (1.3 ± 0.3) X 10-1^ MeV-l (47) 

(Here and below the errors are only experimental, and do not include the uncertainty due 
to higher-order corrections.) For the //-independent ratio this gives: 

4f = (-1.5 ± 0.3) X 10-1° MeV-\ (48) 
^0 ° 
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At this value of p the Coulomb parameter t] = 0.043. The correction proportional to t] in 
Eq. (39) also includes a factor of l/{ap) —0.13. The Coulomb correction is only 3 percent, 
smaller than the uncertainties in the measurement and higher-order effects in EFT(;<^). 

Equation (48) may be used to predict the scattering asymmetry at the higher energy of 



45 MeV, yielding: 



Af{E = 45 MeV) = (-1.7 ± 0.4) x 10" 



(49) 



The two extant low-energy data are thus consistent with a leading-order EFT(^) analysis 
within their combined uncertainties. This is really nothing more than the statement that 
at these energies the asymmetry is scaling with the center-of-mass momentum — as already 
observed in Ref. |26j . 

It should, however, be noted that the center-of-mass momentum for the second experiment 
is already larger than m^r. Sub-leading corrections could therefore be large. A crude estimate 
of these effects can be obt ained by using Eq. (23), which includes the effects proportional to 



(but see also Sec. 



Ill C ) due to strong rescattering. This yields (with 



2.73 fm): 



-1.1 ±0.3) X 10"^" MeV 



(50) 



The shift of ~ 30% with respect to the leading-order value (48) is entirely consistent with 



the expansion parameter of EFT(/). The prediction for the higher-energy datum is now 



Af{E = 45MeV) = (-2.6 ± 0.6) x 10" 



(51) 



In this case the shift is more than 50% of the leading-order value (49), suggesting that the 



point at 45 MeV is indeed too high for profitable application of EFT(/'). The large (par- 
tial) 0{Q) correction computed here suggests that we can anticipate significant additional 
corrections at next-to-next-to-leading order. Given the presence of these corrections, as well 



as the experimental error, there is no real tension between (51) and (46). The large NLO 



correction would, though, seem to imply that the agreement between the datum of Ref. 



and the LO prediction (49) is fortuitous. 



Measurements of the neutron's spin rotation as it passes through parahydrogen have been 
proposed, e.g., in Ref. [M]. The hope here is to extract the longitudinal analyzing power 
of n + p scattering. The thermal energies at which these experiments take place are ideal 
for EFT(/). The leading-order pionless EFT prediction is given in Eq. (24) and depends 



upon the coefficients C^^^'^ ^(A/=of°'*' ^Iai=2)°'' ^ ^''''^^ once low-energy data 

is available constraints on EFT(/) coefficients will result. 



V. CONCLUSION AND OUTLOOK 

We have presented the leading-order low-energy prediction for nn, pp, and np longitudinal 
asymmetries. They depend on five different parameters, but one asymmetry measurement 
each in n + n and p + p, as well as two at different energies in n -|- p, would allow the 
extraction of four of the five parameters. We determined, in agreement with the findings of 
earlier authors, that the Coulomb corrections to pp scattering are not significant at leading 
order for the energies at which these measurements are made. Finally, we showed that when 
the effective range is taken to scale as Q^, the running of the leading order weak interaction 
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coefficients mimics that of Co of the strong interaction, while the running of the next-to- 
leading-order coefficients is expected to mimic that of C2 of the strong interaction. This is a 
simple consequence of the mixing of the S'-wave side of the parity violating operators with 
the bubble-sum enhancement of strong S'-wave scattering. 

Our EFT(/') calculations presented here, which implement a systematic power counting 
scheme, show the consistency of the Danilov hypothesis that (at least for energies < m^/M) 
the dominant energy dependence in parity-violating A^A^ observables arises from the large 
NN scattering lengths in the parity-conserving sector. Effects of energy (or momentum) 
dependence in the parity-violating A^A^ interaction, as well as those due to PV mixing be- 
tween other partial waves, constitute higher-order effects in EFT(;<') which are accompanied 
by additional unknown parameters. 

At leading order in EFT(;f) there are only five independent PV A^A^ operators [HU] . 
This means that five independent measurements will serve to pin down the leading-order 
Lagrangian. Equation ^ is one way to write the five terms of the LO PV Lagrangian in 
EFT(/). It is equivalent to the previously published form ( [5| , with the matching computed 
in detail in the subsequent Appendix. The Lagrangian ^ has the advantage of being 
written in an operator basis where each coefficient contributes to one and only one partial- 
wave transition. ^ 

For instance, the longitudinal asymmetry in pp scattering probes the coefficients 

^(zf/=if°^' ^^"^ ^\li=i)^^ ■ formulae (13) and (23) encode the specific combination in 

which the coefficients associated with different isospm transitions appear. Existing experi- 
mental data can be used to extract this combination — admittedly with large error bars. A 
lower-energy pp experiment with high precision would be a useful development. 

Meanwhile the partial-wave transitions C^'^^-'p^) ^ ctAi=Q)°\ (^(Jrif"^' C^'^^-'^^^ are 
probed in the np longitudinal asymmetry. In principle a detailed study of the energy de- 
pendence of this asymmetry could allow the extraction of C(^'S'i-^^i) and the particular linear 
combination of the other three coefficients relevant for np scattering. 

Finally, the LO EFT(7<') prediction for the longitudinal analyzing power of n + n scatter- 
ing depends upon the coefficients C^^liZl)''\ and Cf^^^Jaf^^ but in a different linear 
combination to that appearing in the prediction for . Given the difficulties inherent in 
such an experiment it seems more productive to focus on the asymmetry 'm.n + d scattering 
at low energies (see also Ref. [32]) and perform the necessary three-body calculations for 
the interpretation of that asymmetry (see Ref. [65t [66] for examples in the parity-conserving 
sector) within the consistent EFT(/) framework for parity- violating A^A^ scattering laid out 
here. 

In any calculation of PV A^A^ observables it is important to treat the PV and PC A^A'^ 
interactions consistently. In particular, care must be taken that operators and wave functions 
used in the same calculation are evaluated using compatible schemes and subtraction points. 
Use of the AV18 A^A^ potential to evaluate matrix elements of the short-range operators in 
Eq. (g |30l [32] represents a significant mismatch in this regard, and cannot be considered 
a systematic EFT(;<^) calculation. 

In Sec. Ill C we emphasized the importance of maintaining a consistent power counting for 
both the weak and strong parts of the Lagrangian. A consistent calculation in EFT(/')can be 
carried out assuming either that r scales as — in which case range corrections are treated 
perturbatively — or that it scales as in which case a dibaryon formalism is necessary. 
The most appropriate choice should be revealed by seeing which (possibly higher-order) 
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predictions provide a better explanation of the data. 

The PV coefficients for NN asymmetries are presently experimentally underconstrained, 
so it is necessary to use electromagnetic reactions in the NN system to probe additional 
linear combinations of the five PV parameters. At lowest order in EFT(/), both np — >■ d^y 
and the anapole moment of the deuteron depend only upon a single coefficient (C^^'^^"^^^-*), 
and so serve to disentangle this coefficient from the linear combination involved in other 
np processes. These have been computed in Ref. [67]. Measurements of the asymmetry 
in np —>■ dj are presently consistent with zero [371 EH], but improvements by an order of 
magnitude are expected [3HI EH] • 

A further constraint on the five PV parameters is potentially available from circularly 
polarized photon-deuteron breakup (or the inverse reaction). Experimentally, results are 
presently consistent with zero [ID]. The development of high intensity free electron lasers 
to produce circularly polarized photons has led to proposals (e.g., Ref. [391 EO]) to perform 
this measurement if the necessary luminosity can be achieved. The PV parameters involved 
in the LO EFT(^)prediction are C^'^^-'^^), C^^fjZo^'\ and Cf^^^Jsf ^ This system has been 
discussed in Ref. [7T]. A partial LO calculation has recently been reported [72] . 

We have presented a model-independent set of operators and coefficients with which 
low-energy PV observables can be described and compared, emphasizing the utility of the 
partial-wave basis. Such a treatment conveys significant advantages in our efforts to under- 
stand manifestations of parity violation in few-nucleon systems. Further calculations and 
experiments which use the EFT(/) framework to map out the landscape of possible exper- 
iments in two-, three-, and four-body systems that are pertinent to parity violation would 
be very useful. 
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VI. APPENDIX 

In this appendix we discuss the matching between the Weinberg basis Lagrange density 
of Eq. ([T]) and the partial wave basis Lagrange density of Eq. (|2|. One method for matching 
the coefficients in these two different bases is to use Fierz rearrangement identities. Another 
method is to use the Lorentz structures with their nucleon spin and isospin indices explicit, 
and employ orthogonality and completeness of the operators in order to isolate one set of 
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operators in terms of the other. The latter is easy to implement using a Mathematica [73] 
code with the HighEnergyPhysics 'FeynCalc' package [71] to perform SU(2) manipulations. 

As a simple example we consider the lowest order terms of the strong interaction Lagrange 
density for the two nucleon system. A useful Fierz identity is 

where e = ia2 is the totally anti-symmetric Levi-Civita tensor, fi is summed over (/i = 
0, 1, 2, 3), cr^ = (1, a) and i,j, k, and I are summed over spin indices. An identical equation 
serves just as well for the r matrices and the isospin indices (we will use a, b, c, and d for 
these). Putting this together: 

{(J^)ij{(Tfj.)kl{T'')ab{Tu)cd = {'i^ik^jl){'i^ac^bd) = 4:{<^2)ik{<^2) jl{T2) ac{T2)bd ■ 

To obtain an equality involving the Cs operator, (A^"''A^)(A^"''A^), we want delta func- 
tions rather than (T2's and T2^s on the right hand side, so act on both sides with 
(o"2)i'i(o-2)n'(7-2)a'a(i"2)dd' and obtalu 

{<y2<y^)i'j{<ytJ.(^2)kl'{T2r'')a'bir,,T2)cd' = 4:Si>kSjl>6a'cSbd' ■ (52) 

Contracting this with Nl^Ni'a'Nj,^,Njb - the nucleon operators with their spin and isospin 
indices explicit - and noticing that diagonal terms with N'^a2T2N and N'^a2(yiT2TjN are 
disallowed by the Pauli principle, we obtain the decomposition of the operator associated 
with C5 in terms of the operators in the partial wave basis: 

(Art AT) (Art AT) = -2{N^PXSo)Ny{N'^PXSo)N) - 2{N^ PiCSi)Ny {N'' P,CSi)N) , 

where a and i are now summed over (1,2,3) as in the projection operators introduced in 
section |Tll 

To obtain the form {N^ aiN){N'^ aiN) , associated with Ct, out of the Fierz identity, we 
need to appropriately insert not only a2 and T2 (to get to delta functions) but cr^i's on the 
right-hand side as well. Contracting 



nd' 



on both sides of Eq. (52) yields 



A{N^a^N){N^aiN) = {N*a^a2r2N^){N^cr2cr,T2N) - {N*a,a2T2rjN^){N^(r2crirjr2N) - 

{N*aia2(yjT2N^){N'^aja2cri T2N) + {N*a,a2ajT2TkN^){N''aja20^TkT2N).{hZ) 

But this can be considerably simplified because 

{N''a,a2T2T,N^){N^(J2(T^T,T2N) = ; 

{N*aia2(Tjr2N^){N^aja2aiT2N) = 2{N*aia2T2N^)iN'^a2aiT2N) ; 
{N*a.,a2a,T2TkN^){N^cr,a2a,rkT2N) = 3{N^a2nT2N^){Na2T2nN), 

(54) 

so that 

(Art a, AT) (Art a, AT) = -2(Ar^P„(i^o)Ar)t(Ar^P,(i^o)Ar) + 6(Ar^P,(3^i)Ar)t(Ar^P,(3^,)A^), 

(55) 
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FIG. 2: Assignment of spin indices {i,j,k,l), isospin indices {a,b,c,d), and momentum labels for 
purposes of matching the partial wave basis operators to the strong Weinberg and weak Girlanda 
operators. Figure created using JaxoDraw |75| . 



-Cs^ik Sac S'jl ^bd + Cs^jk he ^il ^ad — Cxic^ A)ik ^ac {(^ a) jl ^bd + Ct {<^A)jk he (c"A)i/ ^, 



which completes the decomposition of the Cs and Ct operators in terms of the Cq' and Cq° 
operators. 

An easier procedure for obtaining one set of basis coefficients in terms of the other is to use 
orthogonality and completeness of the operator sets. We will illustrate this using, again, the 
leading-order strong-interaction terms. Making the spin indices and isospin indices 

(a,b...) explicit (referring to Fig. |2]) and including all possible nucleon assignments, 

ad 

^ C^0^°\TAT2)ab i(^2)ij {T2TA)ed {(^2) kl + ^ Cf^'^ (^2)06 {(^B(^2)ij (r2)cd (o-20'i?)feZ , (56) 

where summation over A,B = 1,2,3 is implied. Contracting both sides with the first 
structure, 6ik Sac Sji Sm, yields: 

Cj'^°^ + Cl'^'^ = 2Cs - 2Ct . 

A second equation is found by contracting both sides with the third structure, 
{o-A)ik Sac {(^A)ji Sbd, yielding: 

3C5'^°^ - Cj'^'^ = 2Cs - IOCt ■ 

Solving for the partial wave coefficients yields the relationships given in Section [TTJ 

Now consider the two weak-interaction bases from Eq. ([s]) and Eq. (|6|. Note that the op- 
erators in Eq. ^ are explicitly hermitian, but not symmetric under interchange of outgoing 
(or incoming) particles. On the other hand, the basis used in Eq. ^ is symmetric under 
interchange of particles, but each is not its own hermitian conjugate. This is important to 
remember when comparing coefficients. 

Even without Fierzing, inspection of the operators suggests that not all of the partial 
wave operators are involved in the decomposition of, say, the operator associated with Qi 
in Eq. ([s]). But no orthogonality need be assumed because it is easy to verify. The starting 
point is (momenta from Fig. pi): 



4C('^i"'^^)(a^a2)i,(a2)H(r2)a6(r2)cd(2pA)+4C('^^-'^^)(a2aA)H(a2)ii(r2)a6(r2)ed(2p'^) 

+ 4:C'l^liZQ)°\o-2)ijio-2CrA)kliTBT2)abir2TB)cdi2pA) 
+ 4:C'l^AlZo)°\'^2)kl{o-AO-2)ij{rBr2)ab{T2rB)cd{2pA) 
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+ 4C(Aj=i) {(T2)ij{(^2(TA)kl{TBT2)ab{r2Tc)cde3Bc{-'2iPA) 

+ 4^1/^11^"^ ((T2)M((XA(X2)ij(r2rB)cd(rcr2)a6e3Bc(-2ip'^) 

+ 4C(V=2) i(^2)ij{(^20-A)kl{rBT2)ab{r2Tc)cdlBc{'^PA) 
+ ^C(AlZl)''\(^2)kli(^ACr2)ijir2rB)cdirc'r2)ab'lBci'2p'A) 
+ AC^^^''^^'\aACr2)ijicr2ac)klir2)abir2Ts)cd<^ABc{-'^WB) 

+ 4:C^''^^~''^'^\a20-A)ki{crco-2)ij{r2)cd{r3T2)abeABc{-'^ip'B) 
= 2gi{{aA)ikSji - 5ikiaA)ji)ip + p')ASacSbd 

- '2gi{{aA)jkSii - Sjk{crA)il){p - p')Ahc5ad 

- '2Ql(^ABc{(^A)ik{<yc)jl{w' - ip)B5ac5bd 
+ 2gi€ABc{(^A)jk{crc)il{-¥ - ip)BSbcSad 

- 2Q2eABc{(^A)ik{(yc)ji{w - ip)B{{n)achd + ^ac{'rz)bd) 
+ 2g2eABc{(yA)jk{(yc)a{-w' - ip)B{{rz)bc5ad + hc{rz)ad) 

- 2Q5eABc{(^A)ik{ip - ip)B{(Tc)jl^DE{rD)ac{rE)bd 
+ 2Q^eABc{(yA)jk{-ip - ip) B{(yc)illDE{TD)bc{TE)ad 
+ 2gQ5ik{(yA)jl{-ip' + ip)A^BC1,{'TB)ac{rc)bd 

- 20(^5 jk{(y A)ii{ip +ip)AeBC?.{TB)bc{rc)ad ■ (57) 

Again, are the spin indices and a,b,... the isospin indices. The factor of 4 on 

the left hand side comes from all possible nucleon assignments. The coefficients are real. 
Contracting both sides with one operator (in, for example, the partial wave basis) at a 
time yields a set of equations involving coefficients only. For example, contracting with 
(crAO-2)ij(cr2)w(T"2)af,(^2)cd(2p^) yiclds C^^^^"'^!) = liGi- Gi). Provided both sets of opera- 
tors are minimal and complete (at the order desired), there will be a unique solution. The 
procedure is systematized using Mathematica [73] and the results are provided in Section [IT} 
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